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Abstract 

A model system of equations has been derived to describe a toroidally 
rotating tokamak plasma, unstable to Resistive Wall Modes (RWMs) and 
metastable to Neoclassical Tearing Modes (NTMs), using a linear RWM 
model and a nonlinear NTM model. If no wall is present, the NTM 
growth shows the typical threshold/saturation island widths, whereas a 
linearly unstable kink mode grows exponentially in this model plasma sys- 
tem. When a resistive wall is present, the growth of the linearly unstable 
RWM is accelerated by an unstable island: a form of coupled RWM-NTM 
mode. Crucially, this coupled system has no threshold island width, giv- 
ing the impression of a triggerless NTM, observed in high beta tokamak 
discharges. In addition, increasing plasma rotation at the island location 
can mitigate its growth, but does not restore the threshold width. 



1 Introduction 

To operate ITER [T] in Advanced Tokamak (AT) scenarios, which are desir- 
able for high fusion gain, control of performance limiting magnetohydrodynamic 
(MHD) instabilities must be achieved. AT scenarios are characterised by high 
bootstrap fraction and high 13^ = /3{%)/{I /aB), /3 = 2fio{p) /{B'^), (where (p) 
is volume averaged pressure, B{T) is external magnetic field, /(MA) is toroidal 
current, a(m) is minor radius) 2 . Neoclassical Tearing Modes (NTMs) are one 
such performance limiting instability 3^ , with a detrimental effect on the achiev- 
able Pn- NTMs are unstable if two criteria are satisfied: (3 must be sufficiently 
high, and the plasma must typically be "seeded" with a magnetic island pro- 
duced by another activity |4, ,5, , such as sawteeth [SI [71 [S] , fishbones [8 or ELMs 
[5]. This initial seed island will only grow, and subsequently saturate, if it ex- 
ceeds a certain critical width [TU]. The flattening of the pressure profile caused 



by the presence of the island, and the consequential removal of the pressure 
gradient-driven bootstrap current in the vicinity of the magnetic island 0-point 
enhances filamentation of the current density, providing a feedback mechanism 
for continued island growth. However, NTMs which do not appear to have grown 
from a seed island, so-called triggerless NTMs, have also been observed exper- 
imentally [HI im HH IS], usually associated with high /?. Another instability 
which is responsible for limiting performance at Pn > ^^O"™"" is the Resis- 
tive Wall Mode (RWM); a branch of the long- wavelength external kink mode, 
driven in a tokamak by pressure but which has a growth rate slowed by the 
presence of finite resistivity in surrounding conducting walls. RWMs are slowly 
growing on the timescale of the wall diffusion time, and have a complex rela- 
tionship with plasma rotation. Analytical models predict that a modest amount 
of plasma rotation is sufficient to stabilise RWM growth [131 EH EH IH] ; and 
both experiment and theory show that RWM drag can damp plasma rotation 
[131 HZ] ■ Here, a model system of equations is derived containing linear RWM 
and nonlinear NTM models and extended to include self-consistent nonlinear 
interactions with toroidal plasma rotation. These equations are examined in a 
limiting case relevant to the tokamak and applied as a possible model for explor- 
ing the phenomenon of triggerless NTMs. Section[2]will outline the derivation of 
the model system considered. This is followed by an examination of the stability 
of the RWM-NTM coupled mode and its dependence upon model parameters 
in Section [3l conclusions are discussed in Section [H 



2 Theoretical Model 

Suppose an ideal, toroidally rotating plasma is surrounded by a thin wall at 
r = with finite conductivity a^, and contains a resistive internal rational 
layer at radius r — Vs where a magnetic island can form. The outer ideal 
regions of the plasma are treated linearly, whereas nonlinear effects associated 
with magnetic islands will be retained in the resistive plasma of the layer. 

Due to skin currents in the resistive wall and at the rational surface, the 
parallel component of the vector potential has a discontinuous radial derivative 
at these locations. These discontinuities are parametrised by and re- 
spectively. In a linear cylindrical model employing a complex representation of 
the perturbed fields, they are related through an expression of the form |18j 

where e and 5 are related to the stability properties of the plasma in the limit 
of no wall and a superconducting wall respectively. The parameters u, are 
derived from ideal MHD outside the wall and rational layer. They must be 
matched to the quantities derived from the current perturbations flowing in the 
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layer/wall and their effect upon the magnetic flux: 



^L,w — w— 

ipc or 



(2) 



Here is the complex representation of the magnetic flux function derived 
from linear ideal MHD. When no wall is present, A^, = and for an ideally 
unstable plasma, e is positive and small. For an ideal superconducting wall 
at r = r^, A„ — ^ oo, and (5 > provides stability to the kink mode. Thus 
small e, 5 > describe a plasma susceptible to a RWM. Whilst we solve linear 
equations in the outer region, we must anticipate a need to match to nonlinear 
layer solutions at the rational surface. Nevertheless, we assume we can separate 
the time dependence and write the outer region solution in the form: 

V'c(r,t) =Vi(r)e'«e-^P''* (3) 

where p{t) = "f{t) — iu}(t), with 7 the instantaneous growth rate, uj the mode 
frequency relative to the wall and ^ = m9 — n(p the helical angle defined in terms 
of the poloidal {9) and toroidal {(j>) angles. It will be assumed that the toroidal 
mode number n = 1, with m the poloidal mode number. Using Ampere's Law, 
the wall response can be written as 

= pr (4) 

where t = r^/Tr and — iioda^rm is the wall diffusion time, with the wall 
thickness d « r^, the minor radius at which the wall is located. The instanta- 
neous complex growth rate p is normalised to = crfiQars the resistive plasma 
diffusion time where a is the plasma conductivity at the rational surface. Com- 
bining these results relates the layer response A^ to the wall properties and the 
instantaneous complex growth rate |18| : 

A. ^ ^ (5) 

-e + TP 

In the ideal linear region, a complex representation for the magnetic flux 
ensures that relation ([T]) holds. The physical flux is denoted ipL = Relipc]- In 
the rational layer, we work with this physical flux and assume a magnetic island 
exists at the rational surface. Let us transform into the frame of reference where 
the rational surface is at rest, so that the wall rotates. Analogous to the wall. 
Ampere's law describes the current perturbations 

-Q^ = J|| (6) 

where J|| is the component of current perturbation parallel to the magnetic 
field. There are two key contributions to J|| to consider: an inductive component 
proportional to dip L/dt and the perturbation in the bootstrap current J^s caused 
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by the pressure flattening inside the island. We assume that these current 
perturbations are locahsed within the layer and integrate across that layer to 
derive the model equation: 



dr 



2afj,ow^^^ - 2fj,awJbs (7) 



where w is the half-width of the magnetic island and a is the plasma conduc- 
tivity. The discontinuity on the left of Equation ([t]), caused by the current 
sheet in the layer, can be expressed in terms of as follows. Let V'l be of 
the form ipL = (pcos{a + C), where d = ^ — J ujdt, C is a phase factor, and 
uj ^ u ~ VLl, such that is the toroidal plasma rotation frequency of the 
plasma at r = rs with respect to a stationary wall. Matching the ideal solution 
provides Rel^tle^ '^'^^ = (pcos{C) and Im[ip]e-I' = (psin{() (we assume that the 
toroidal mode number n=l and ip is the amplitude of the magnetic flux function 
in Equation ([s])). Thus, combining the results yields 

((^/r,)(i?e[Ai] cos(a + C) - /m[Ai] sin(a + C)) = 

2afiow{^ cos(a + C) + t^(psin(a + C)) - 2fiowJi,s (8) 



Multiplying by cos(q; -|- C) or sin(Q; + () and then integrating over a yields 

dip 
dt 



equations for uj and The island half width w is related to (f by 



where Lg is the shear length scale, q the safety factor, R the major radius and 
Bg the poloidal magnetic field. Normalising the island half width w to the minor 
radius, and time to we obtain 

4w^Re[AL] + -(l-^) (10) 



w 

2w{uj-nL)+Im[AL]=0 (11) 

where w refers to time derivatives with respect tot = t/rr and Jbs — '^^^^^ cos(a-|- 
^), where f3g is the poloidal beta, e the inverse aspect ratio, Lp the pres- 
sure length scale, and Wc the seed island threshold width. We have defined 
/3 = ^^^l^' , s the magnetic shear. Note that 7 = 2w/w, uj and fi/, are nor- 
malised to Tr- We have also introduced a heuristic threshold factor (1 — w^/w^), 
which could be attributed to the polarisation current effects in the plasma jl9] . 
We have neglected geometrical factors to yield a simpler model which never- 



theless retains the essential physics. Equation ( 10 ) is the same relation as that 
found in |20| when rotation is not considered. 

To close the system, we require a torque balance relation. Ideal plasma is 
torque- free [21]; the torque exerted on the plasma is a delta function at the 
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rational surface S{r = rg). 

Consider the perturbed MHD momentum conservation equation in the non- 
ideal layer, averaged over the flux surface. 



{6J X 6B) + p^V^v = 



(12) 



where angled brackets denote the average. We assume that the processes under 
consideration occur over many diffusion times, and pressure is constant across 
the thin layer. Here p is the plasma density and /i is the plasma viscosity. The 
pressure gradient is neglected as it will not contribute once we integrate across 
the layer (as follows). Assume that v is continuous across the layer but Vw 
possesses a discontinuity at r = , due to the localised torque at that location. 
For a perturbation (5B, V x (5B = — V^f/'Lb. Only the toroidal component is 



required. Integrating (12) over the rational layer 



-smla 



C) 



SipL 



pp. 



Sr 



(13) 



The discontinuity in the radial derivative of ipL = Rs['<Pc] is evaluated as above. 
Integrating over a, the terms in i?e[AL] disappear. Let = flR and normalise 
as before to derive 



sn 

Sr 



A 



(14) 



such that A = s^e^a^TvTr/512rgq^T^,with ry the momentum confinement time, 
a the minor radius, g=2 the safety factor at r = r^, the Alfven time. The 
electromagnetic torque is finite when /m[Ai] 0, and always acts to damp the 
plasma rotation [22]. Im[Ai,] is related to the phase factor ( [18]: the torque 

only acts on the plasma in the rational layer when ^0. 

In the outer ideal plasma regions, the momentum equation is simply 







(15) 



The linear rotation profile in the ideal plasma regions, r < and < r < a, is 
continuous &i r = but has a discontinuity in the first radial derivative there, 
caused by the torque. Imposing no slip boundary conditions at the outer edge 
of the plasma r = a (viscous drag effects) and a driving force at the inner edge 



at r = that maintains 



dn I 



-Qo/a, the toroidal rotation frequency profile 



can be constructed in the region of ideal plasma. 

ni(r) = A(a-r,) + ^(r, -r) < r < 
^2(1^) = X{a — r); < r < a 



(16) 
(17) 



for constant A > 0. The discontinuity in the radial derivative is calculated to be 
—A + rJo/o- Hence, at r = Ts, using Equation (14), the damping of the toroidal 
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plasma rotation frequency is found to be 



(18) 



Numerical Solutions For Coupled RWM-NTM 
Mode Stability 



The system of nonlinear equations ([5|), ( 10 1, (111, and ( 18 1 can be solved numer- 
ically and integrated, with r w 10 (a — 1.0m, rg = 0.5m, = lkeV,lnA = 
20). If there is no wall surrounding the plasma, then r = 0. In this situation, 
the RWM is not a solution (an infinite growth rate is obtained, which can be 
interpreted as the ideal kink mode). However, the NTM is a solution in this 
limit (Figure [ij. A sufficiently high /3 must be achieved for the NTM to become 
unstable, as well as a seed island with a width w > Wc- An NTM will saturate 
at a large island width, proving detrimental to plasma confinement. 
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Figure 1: No RWM exists when the resistive wall is removed: only the NTM is 
a solution. The dependence of the island growth rate on /3 and Wc is shown. 

When r is finite, the wall is considered to be either close-fitting, thin, or 
with a low conductivity. Both RWM and NTM branches are present. 

In the case with zero /3, the NTM has no drive and we have a conventional 
RWM: wall- locked RWM acting to slow plasma rotation which also eventually 
locks to the wall. The mode grows exponentially from any seed width since the 



6 



RWM dispersion relation ([ij is linear. Increasing /3 destabilises the NTM for 
sufficiently large values, and the RWM takes on a different character. We now 
examine the situation where a seed island has a width smaller than the NTM 
threshold width Wc, such that the NTM solution is stable. We observe that a 
seed island of any size is able to grow, as in Figure [2^) . In Figure [2^) , the seed 
island grows steadily until it reaches width Wc- At this point, the island growth 
rate increases dramatically for high /3. If f3 is reduced, then the NTM drive 
is reduced and the increase in island growth is either reduced or insignificant. 
The plasma rotation frequency at r = is damped by torque exerted by the 
unstable RWM, until the island width reaches Wc (Figure [2];)). The plasma 
then briefly spins up, over a time scale of 3-4 ms, then slows to lock to the wall. 
The mode locks to the plasma simultaneously with the plasma spin up and the 
island crossing the NTM threshold width. In this case, increasing f3 increases 
the size of the spin-up. 

The island evolution can also be influenced by the amount of momentum 
injected at the inner edge of the plasma, as characterised by f2o- This is shown 
in Figure |3] Increasing the initial plasma rotation freqency will damp the island 
growth rate when the island width become larger than Wc, although the island 
still grows below that threshold. Increasing the initial plasma rotation frequency 
prevents the plasma spin up and also prevents the mode locking to the plasma. 
For sufficient initial plasma rotation, the mode remains locked to the wall. 

4 Conclusion 

The features described above of a mode with a strong /? dependence and the 
possibility for an island to grow from any seed width, with an increase in growth 
rate after the island achieves the threshold width, are features associated with 
the triggerless NTM phenomenon. 

Although strictly we observe this coupled RWM-NTM mode, it would appear 
to an experimentalist as a triggerless NTM. It can also be seen that increasing 
the initial plasma rotation frequency can mitigate the island growth: this could 
be due to the mitigating effect of plasma rotation on the RWM growth. It should 
be noted that increasing plasma rotation frequency does not prevent very small 
seed islands from growing, but only reduces their growth rate. The plasma spin 
up observed as the island width approaches the threshold for low torque input 
could be used to test this RWM-NTM coupled mode theory. 

The phenomenon of triggerless NTMs near j5'}^~^°-^\ the ideal no- wall /3 
limit, can be explained by invoking a pole in A (which here is equal to Ai/r^). 
As /37V — >■ A — >• oo, destabilising a classical tearing mode, which 

in turn destabilises the neoclassical tearing mode when at a sufficiently large 
amplitude [T2]. In our model, this would be equivalent to letting A^u and 
e — )• 0, corresponding to the no-wall ideal MHD stability boundary. In the model 
described here, the RWM is destabihsed for Pn > /3]^°-'"'^" but the pole in A^ 
is resolved by dependence on the wall response. The physics interpretation is 
therefore different to that provided in [T^] . The features of the triggerless NTM 
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appear in this system, although it is the RWM which exhibits them; the true 
NTM branch retains its threshold and saturation widths. Finally we note that 
our RWM branch does not saturate ar large island widths because we have not 
retained the nonlinear RWM physics. 
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Figure 2: Here w{t = 0) = 0.005 which is considerably smaher than the critical 
seed island width — 0.02. The dependence of the island growth, mode 
frequency and plasma rotation frequency at the rational surface on /3 is shown: 
increasing (3 (which increases the NTIvr drive) allows the island to grow faster 
and to a larger size. The rotation spin-up is however reduced and eventually 
removed. 
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Figure 3: w^t = 0) = 0.005, /3 = 1.0 a) demonstrates that the initial plasma 
rotation frequency has a noticeable effect on the island evolution, whilst b), 
c) and d) show the mode and rotation frequencies for differing values of fio- 
rio = 10000 corresponds to the plasma rotating at about 5% of the plasma 
sound speed. 
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